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1.  INTRODUCTION 


As  part  of  the  fatigue  life  assessment  of  the  R.A.A.F.  Macchi  MB326H  aircraft,  a full 
scale  fatigue  test  was  carried  out  by  the  manufacturer,  Aeronautica  Macchi,  S.P.A.,  Varese 
(Ref.  I).  \\ 

Since  ihenTlesting  techniques  have  been  developed  to  allow  more  realistic  load  sequences. 
To  take  advantage  of  this,  the  design  of  such  test  sequences  must  also  be  improved  and  for  several 
years  large  scale  tests  at  ARL  have  relied  on  loads  obtained  from  flight  histories  transformed 
as  little  as  possible.  In  this  way  the  probability  structure  of  flight  sequences  can  be  restored. 

Until  now  however,  the  design  of  test  histories  has  been  almost  entirely  manual.  This  note 
describes  the  generation  of  histories  by  use  of  automatic  data  processing  techniques  and  linear 
programming. 

The  example  used  in  this  note  is  the  Macchi  centre  section  spar  boom  which  is  being  tested 
under  the  design  load  history  to  supplement  earlier  data  obtained  under  less  realistic  load 
sequencejjsiRef.  2). 

2.  DATA 


2.1  Fatigue  Meter  Spectrum 

The  KM  spectrum  for  this  analysis  was  obtained  from  magnetic  tape  *1  containing  R.A.A.F. 
FM  data.  The  data  was  from  flights  at  No.  2 Flying  Training  School  (2BFTS)  Pearce,  during 
1973,  1974  and  1975,  on  aircraft  fitted  with  70  gallon  fuel  tanks  and  no  stores. 

It  comprised  15,748  flights  made  up  of  9 types  of  flying  (TOF).  Of  these,  2 were  flown  so 
infrequently  (9  times)  that  they  could  be  safely  ignored.  The  remaining  7 were  analysed  (Table  I) 
and  then  scaled  to  200  flights  (Table  2). 


2.2  Flight  Loads 

Loads  were  obtained  from  histories  of  flight  trials  conducted  in  1974  at  Laverton  and 
Williamtown  (Ref.  3,  Table  3). 

The  test  aircraft  was  Macchi  A7-005.  Forty-nine  channels  of  analogue  data  were  each 
sampled  60  times  per  second  by  a multiplexer,  converted  to  digital  form,  and  written  onto 
magnetic  tape  on  a recorder  mounted  in  the  rear  cockpit. 

The  following  variables  were  recorded:  normal,  lateral  and  longitudinal  acceleration: 
roll,  pitch,  and  yaw  rates:  rudder,  aileron,  and  elevator  position;  airspeed;  altitude;  sideslip; 
incidence  and  rudder  force.  The  remaining  channels  were  connected  to  electrical  resistance 
strain  gauges  on  the  lin  and  the  centre  section  main  spar. 

3.  METHOD 

3. 1 Sequence  of  Flights 

The  method  followed  for  producing  the  flight  sequence  was  that  employed  by  Howard 
(Ref.  4). 

Two  event  frequency  tables  were  produced  by  analysis  of  the  Pearce  flight  data.  The  first 
(Table  4)  shows  the  frequency  with  which  the  subject  mission  is  followed  by  a different  mission 
(transition  frequency).  The  second  (Table  5)  shows  the  frequency  of  occurrence  of  various 
run  lengths  (number  of  times  the  mission  is  flown  before  a different  mission  occurs). 

The  ordering  observed  is  similar  to  that  found  for  the  Mirage  1 1 10  by  Howard.  It  will  be 


*1  provided  by  AIRFNG5A I,  DEFAIR. 


noted  that  the  row  sums  m Table  4 must  sometimes  dill'er  from  the  corresponding  column  sums. 
This  follows  from  the  fact  that  only  transitions  to  unlike  missions  are  included  and  that  row 
sums  are  following  missions  whilst  column  sums  count  missions  before  a given  type.  Further 
analysis  indicates  that,  save  for  inevitable  miscounting,  corresponding  row  and  column  sums 
are  mostly  equal.  The  two  exceptions  are  the  counts  for  first  and  last  flight  types  which  may 
differ  by  one. 

The  KM  data  analysed  related  to  1 5,748  flights,  and  had  to  be  scaled  to  200  flights  for  the 
test.  This  was  accomplished  by  multiplying  all  values  in  the  frequency  tables  by  200  I5.74X. 
Fractional  values  in  the  transition  frequency  matrix  were  eliminated  by  rounding  off  to  the 
nearest  integer  The  run  length  table  required  slightly  more  complicated  treatment.  Flights 
representing  fractional  frequencies  of  run  length  were  accumulated  until  an  integral  number  of 
shorter  runs  could  be  formed.  For  example  a frequency  of  0-5  for  a run  length  of  12  and  0-55 
for  a run  length  of  1 1 became  one  run  of  1 1 with  one  (light  carried  to  shorter  run  lengths. 

On  rounding  off  the  transition  frequency  matrix,  a number  of  (lights  were  lost  and  had  to 
be  restored.  The  method  used  is  shown  in  Appendix  I.  The  resultant  matrix  is  shown  on  Tables  6 
and  7.  In  Table  6 the  scaling  and  rounding  has  also  removed  the  differences  between  row  and 
column  sums  [his  situation,  which  is  physically  impossible,  leads  to  the  end  effect  described 
below. 

The  actual  (light  sequence  was  created  using  the  matrices  in  Tables  6 and  7.  This  was  done 
in  the  following  way: 

A TOF  is  selected,  in  this  case  the  least  severe.  TOF  0 (to  ensure  that  the  sequence  does  not 
commence  with  a high  load  application).  A subsequent  TOF  is  then  selected  at  random  from 
the  transition  frequency  table,  the  next  step  is  to  look  at  the  run  length  table  and  choose  a run 
length  at  random  of  this  subsequent  TOF.  Then  back  to  the  sequence  table  to  choose  the  next 
following  TOF  and  so  on.  each  time  subtracting  one  from  the  appropriate  position  in  each  table. 
This  is  repeated  until  the  tables  contain  all  zeros. 

Inherent  in  this  method,  however,  is  an  end  effect.  This  adds  an  extra  (light  into  the  matrix 
at  the  start.  Before  the  tables  are  completely  exhausted  we  find  ourselves  in  a row  containing 
all  zeros  and  therefore  unable  to  select  the  next  TOF  or  proceed.  To  overcome  this  difficulty 
we  backstep  to  the  last  transition  and  subtract  one  from  the  matrix  at  the  point  that  led  into 
the  all  zero  line.  Then  go  back  to  the  beginning  and  start  again.  Though  apparently  fairly  tedious, 
use  of  the  digital  computer  greatly  facilitates  the  effort. 


.4.2  Compression  of  Flight  Trials  Data 

As  fatigue  testing  only  requires  turning  points,  all  other  data  can  be  filtered  out.  The 
compression  of  the  flight  trials  data  was  accomplished  in  two  stages.  The  first  stage  was  to 
extract  data  that  may  be  used  in  a later  fatigue  analysis  of  the  structure,  and  the  second  extraction 
extracts  from  these  data  only  what  is  required  to  provide  load  sequence  in  the  spar  boom  for 
the  fatigue  test. 

STAGE  I 

Fifteen  channels  were  selected  for  data  extraction  (Table  X).  These  channels  have  been 
found  to  be  the  most  important  from  previous  analyses  (Ref.  5 and  6).  Before  a turning  point 
can  be  extracted,  discriminators  need  to  be  chosen  in  order  to  eliminate  negligible  load  changes 
(effectively  the  smallest  range  pairs).  The  discriminator  is  the  required  difference  that  must  be 
exceeded  before  a turning  point  is  registered.  For  example  (see  Fig.  I)  if./  discriminator  and 
among  the  four  turning  points: 

e f g h 
and 

/ h anil  g e 

then  if  g f (I.  g and  /Will  be  excluded  from  the  data. 

■> 


The  discriminator  chosen  were: 


For  normal  acceleration  (/»..): 

For  strain  gauges  on  centre  section: 
(strain  is  linearis  related  to  «..) 

For  strain  gauges  in  the  tin: 


0 25  g 

strain  corresponding  to 
0-25  g 

an  arbitrary  value  of 
50  microstrain 


e 


FIG.  1.  USE  OF  DISCRIMINATOR 


I he  fifteen  channels  selected  were  scanned  by  a computer  program  (MACFXT.F4).  and 
when  a turning  point  occurred  on  any  one  of  the  fifteen  channels  then  the  data  from  all  49  channels 
were  converted  to  engineering  units  (Ref.  7),  scaled  and  written  onto  magnetic  tape.  The  com- 
pessed  data  occupy  three  tapes  entitled  "Macchi  MB  7 26 FI  Flight  Trials  Stage  3 Edited  Data  , 
I ngineering  Units”  tapes  I.  2,  V 


STAGE  2 

t rom  each  ol  the  I diled  flight  tapes,  the  Normal  Acceleration  n:  was  extracted  and  passed 
through  .i  turning  point  routine  t If  RNPI  I 4)  to  produce  nine  tiles  containing  turning  points 
of  ii;  each  representing  a different  history.  Due  to  the  small  number  of  histories  available,  a 
number  of  quasi-flights  was  generated  by  scaling  every  turning  point  on  the  nine  flights  by  the 
relation 

A ( I)  I (Ref.  4) 

where  «„  is  the  scaled  I oad  f actor,  n is  the  original  Load  f actor  and  A has  the  value  0-8  to 
I 4 incremented  by  0 I 

Ihe  values  ol  A selected  have  been  used  previously  (Ref.  4)  and  were  found  to  be  suitable 
in  this  project  Ihe  scaling  produces  the  same  type  of  flying  but  with  allowances  for  different 
pilots  flying  the  same  mission 

I his  produced  7 flights  for  each  of  the  9 flight  histories,  giving  a total  of  63  flights  (Fig.  2). 

A Macchi  MB  32611  fatigue  meter  tvpc  count  was  done  on  the  63  flights  to  produce  an 
FM  spectrum  for  each  t fable  9) 


5.0 

0 


Number  of  turning  points 


FIG.  2.  THREE  TYPICAL  FLIGHTS 


3.3  Creation  of  a Test  f light  Sequence  Tape 

Analysis  of  the  fatigue  meter  data  produced  a series  of  200  (light  numbers  in  which  7 TOF's 
occured  in  a sequence  typical  of  squadron  usage.  Nine  (light  trials  (lights  were  processed  as 
above  into  63  histories.  To  create  the  test  (light  sequence  (ape.  a (light  history  had  to  be  matched 
to  each  of  the  200  (light  numbers  This  would  be  a simple  matter  if  each  of  the  (light  histories 
represented  a TO!  classification.  Fhis  is  not  the  case  ( Ref.  6)  and  so  the  TOF's  and  the  9 (light 
trials  flights  were  grouped. 

The  grouping  involved  placing  of  histories  and  TOF 's  of  similar  flight  descriptions  into 
separate  groups.  On  performing  the  linear  programming  to  match  the  histories  with  TOF’s, 
no  feasible  solution  was  possible  until  the  number  of  groups  was  reduced  to  2 (Table  10). 

An  FM  spectrum  wa  • drawn  for  the  2 ty pes  of  fly  ing  ( Tabli  II).  It  was  then  required  to  find 
out  what  combination  of  the  flight  trials  (including  the  quasi-flights)  will  produce  the  same  FM 
spectrum  and  at  the  same  time  maintain  the  number  of  each  type,  as  shown  in  table  II.  This 
was  achieved  by  linear  programming  for  each  type.  (See  Appendix  2 for  worked  example.) 
Flight  I3SI  was  unable  to  be  matched  up  and  was  slightly  altered  to  achieve  a match.  The  new 
(light  is  coded  I38B.  Lhe  results  are  shown  on  Table  12.  As  can  be  seen  from  Tables  1 1 and  12 
a good  match  has  been  achieved.  By  using  the  information  from  Table  10  with  the  results  shown 
in  Table  12,  a computer  program  (MASFQ.F4)  moves  down  the  (light  sequence  list  (Sec.  3.1) 
and,  (light  by  (light,  (using  (light  trials  histories)  writes  onto  magnetic  tape  the  appropriate  history . 

It  is  at  this  point  that  corrections  are  carried  out  (i.e.  inserting  one  6g  level,  subtracting 
two  I -5  g,  adding  one  0-5  g and  subtracting  one  2 • 5 g levels).  A fatigue  meter  type  count 
was  then  done  on  the  complete  tape  of  the  200  (lights  the  results  of  which  appear  on  Table  13 
together  with  the  totals  from  Table  1 1. 

The  final  product  is  a magnetic  tape  containing  200  files,  entitled  “Macchi  fatigue  test 
200  flights  Turning  Points".  Each  file  represents  one  (light.  The  FM  spectrum  of  this  tape 
matches  that  required  from  the  Pearce  I'M  data  (Table  13.  F:ig.  3). 
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Io  achieve  representative  loading,  each  n:  value  on  the  test  tape  must  be  converted  to  load 
units.  I he  effective  critical  stress  sensitivity  lor  the  weakest  section  boom  currently  living  (04a) 
is  92  - 5 MPa  n.  (Re).  X)  I he  end  load  on  the  test  specimen  (an  04  of  smaller  area)  is  determined 
by  increasing  the  stress  m the  ratio  of  boom  areas,  back  calculating  the  expected  strain  at  the 
monitored  strain  gauge  and  then  by  experiment  calculating  the  end  load  against  gauge  response. 
In  symbols,  test  load 

, . -4o4  « critical  t experimental  at  gauge 

/■ill  4|,|  • 

Toi  i <’  gauge  < predicted  at  gauge 

where  A,  a,  < are  boom  areas,  calculated  stress  and  strain  respectively 
The  value  arrived  at  is: 

Load  24-88  n:  Kips 
or  in  SI  units 
Load  1 10-67  n-  A V 

4.  DISCUSSION 

A simulated  sequence  of  200  (lights  representing  a scaled  sequence  of  actual  (lying  en- 
countered under  squadron  usage  has  been  produced.  This  will  be  used  in  the  centre  section  spar 
fatigue  test. 

This  sequence  is  only  one  of  many  possible  combinations.  Other  sequences  could  be  generated 
by  changing  the  starting  point  and  or  using  different  random  data. 

The  technique  employed  requires  very  little  judgement  by  the  operator  and  appears  to  be 
widely  applicable  The  main  drawback  is  that  it  relies  heavily  on  the  existence  of  an  adequate 
set  of  representative  (light  loads  and  fatigue  meter  data,  which  are  not  alway  s available. 

A list  of  computer  programs  used  and  their  functions  can  be  found  in  Appendix  3. 
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Exceedances  per  200  flights 


FIG.  3.  LOAD  SPECTRUM  FOR  BOTH  PEARCE  FATIGUE  METER  DATA  AND  ACTUAL 
FATIGUE  TEST  TAPE. 
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TOF 

1 -5 

Exceedance  Levels  (g) 
0-5  2-5  3-5 

4-5 

6 0 

Number  of 
Missions 

0 

215 

3580 

19983 

12592 

2942 

13 

1992 

5 

820 

12827 

79584 

44054 

8846 

48 

8863 

10 

18 

358 

2918 

1125 

195 

1 

1484 

20 

15 

187 

3317 

1 164 

192 

1 

898 

30 

16 

300 

2219 

1437 

426 

9 

194 

40 

27 

415 

9770 

5681 

953 

6 

2180 

70 

1 19 

714 

3096 

2320 

583 

4 

128 
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TABLE  2 

Above  table  scaled  to  200  flights 


TOF  Exceedance  Levels  (g)  Number  of 


1 -1-5 

-0-5 

2-5 

3-5 

4-5 

60 

Missions 

0 

3 

45 

251 

158 

37 

0 

25 

5 

1 1 

164 

1015 

562 

113 

1 

113 

10 

0 

5 

37 

14 

2 

0 

19 

20 

0 

2 

41 

14 

2 

0 

11 

30 

0 

3 

23 

15 

4 

0 

2 

40 

0 

5 

126 

73 

12 

0 

28 

70 

2 

1 1 

48 

36 

9 

0 

2 

Total  200 


ARDU* 
Recorded 
Hight  No. 

Date 

Pilot 

Location 

1376 

12-2-74 

Dickens 

Laverton 

1377 

12-2-74 

Dickens 

Laverton 

1378 

13-2-74 

Dickens 

Laverton 

1379 

14-2-74 

Dickens 

Laverton 

1381 

15-2-74 

Dickens 

Laverton 

1382 

27-2-74 

Dickens 

Laverton 

1385 

7-3-74 

Doyle 

Laverton 

1386 

8-3-74 

Doyle 

Laverton 

1387 

12-3-74 

Ashbrook 

Laverton 

* Aircraft  Research  and  Development  Unit  — R. A. A. F. 


TABLE  4 

Flight  Transition  Matrix 


TOF 

0 

5 

Following  TOF 

10  20  30 

40 

70 

0 

0 

774 

158 

104 

20 

184 

1 1 

5 

782 

0 

674 

424 

71 

951 

61 

10 

152 

676 

0 

104 

13 

115 

15 

20 

99 

407 

1 10 

0 

8 

96 

3 

30 

15 

81 

14 

7 

0 

17 

2 

40 

197 

973 

100 

78 

21 

0 

12 

70 

5 

51 

18 

8 

3 

19 

0 

TABLE  9 

KM  Spectrum  for  flight  Trials  and  Quasi-Flights 

Macchi  RK  1376  12-2  74 


Hz 

0 8 

Exceedances  for  Each  K Factor 

0-9  10  11  1-2  13 

1 4 

—2-5 

0 

0 

0 

0 

0 

0 

0 

-1-5 

0 

0 

0 

1 

1 

2 

2 

-0-5 

2 

2 

3 

4 

5 

6 

6 

2-5 

12 

14 

15 

15 

15 

16 

16 

3-5 

6 

8 

9 

10 

13 

13 

13 

4-5 

0 

1 

2 

5 
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8 

9 

0 

0 

0 

0 

0 

1 

2 

80 

0 

0 

0 

0 

0 

0 

0 

Macchi  RF  1377  13-2-74 


nz 

0-8 

0-9 

Exceedances  for  Each  K Factor 
10  11  1-2 

1-3 

1 -4 

2-5 

0 

0 

0 

0 

0 

0 

0 

I -5 

0 

0 

0 

1 

3 

3 

3 

0-5 

3 

4 

4 

5 

5 

5 

5 

2-5 

9 

9 

10 

14 

16 

17 

17 

3-5 

7 

9 

II 

11 

11 

11 

11 

4-5 

0 

1 

3 

6 

8 

9 

11 

60 

0 

0 

0 

0 

0 

1 

3 

8 0 

0 

0 

0 

0 

0 

0 

0 

Macchi  Rl  1378  13  2 74 


Hz 

0 8 

0 9 

Exceedances  for  Each  K Factor 
10  ||  12 

1 3 

1 4 

| 

r j 

0 

0 

0 

0 

0 

0 

0 

1 -5 

0 

0 

0 

0 

0 

0 

0 

-0-5 

0 

0 

0 

0 

0 

0 

0 

2-5 

0 

0 

0 

0 

0 

1 

1 

3-5 

0 

0 

0 

0 

0 

0 

0 

4-5 

0 

0 

0 

0 

0 

0 

0 

6 0 

0 

0 

0 

0 

0 

0 

0 

80 

0 

0 

0 

0 

0 

0 

0 

TABLE  9 (continued) 
Macchi  RF  1385  (030)  7-3-74 


Macchi  RF  1386  (31)  8-3-74 


n. 

0 8 

0 9 

Exceedances  for  Each  K Factor 
10  II  1-2 

13 

1 4 

-2-5 

0 

0 

0 

0 

0 

0 

0 

-15 

0 

0 

0 

0 

0 

0 

0 

-0  5 

0 

0 

0 

0 

0 

0 

0 

2-5 

1 

1 

1 

2 

4 

4 

4 

3-5 

0 

0 

0 

0 

1 

1 

1 

4 5 

0 

0 

0 

0 

0 

0 

0 

60 

0 

0 

0 

0 

0 

0 

0 

8 0 

0 

0 

0 

0 

0 

0 

0 

Macch.  RF  1387  (32)  12-3-74 


n. 

0 8 

0 9 

Exceedances  for  Each  K Factor 
10  M 1-2 

1-3 

14 

-2-5 

0 

0 

0 

0 

0 

0 

0 

-1-5 

0 

0 

0 

0 

0 

0 

0 

-0  5 

0 

0 

2 

3 

3 

3 

3 

2-5 

12 

12 

12 

12 

14 

14 

14 

3-5 

6 

7 

10 

12 

15 

15 

15 

4-5 

2 

3 

4 

6 

6 

7 

10 

6 0 

0 

0 

1 

2 

3 

3 

4 

8 0 

0 

0 

0 

0 

0 

0 

1 

Table  9 (continued) 
Macchi  RF  I38B  (26)  15-2-74 


0 8 

Exceedances  for  Each  K Factor 

09  10  II  1-2  1-3 

1-4 

-2-5 

0 

0 

0 

0 

0 

0 

0 

-1-5 

0 

0 

0 

0 

1 

1 

1 

-0  5 

1 

1 

1 

1 

1 

1 

1 

2-5 

II 

II 

19 

19 

19 

19 

19 

3-5 

2 

6 

9 

14 

13 

16 

14 

4-5 

0 

0 

1 

2 

4 

6 

9 

0 

0 

0 

0 

0 

0 

0 

Efl 

0 

0 

0 

0 

0 

0 

0 

TABLE  10 

Flight  Trials  Flights  Corresponding  to 

Pearce  Flights  i 


TOF 

(Pearce) 


0,  5,  30,  70 


Flight  Trials 
Histories 


10,  20,  40 


TABLE  11 

FM  Spectrum  for  the  Two  Types  of  Flying  Done  at  Pearce 


Type 

Exceedance  Levels  (g) 

No.  of 

-1-5 

-0  5 2-5 

3-5 

4-5 

6 0 

flights 

1 

16 

223  1337 

771 

163 

1 

mm 

2 

0 

12  204 

101 

16 

0 

Total 

16 

235  1541 

872 

179 

■ 

200 

APPENDIX  1 
Plight  Sequence  Example 


To  till  in  the  missing  flights  the  following  method  is  used: 

Starting  at  the  bottom  left  hand  corner  of  Table  A I,  compare  row  and  column  totals  (I  and  16) 
with  the  required  totals  from  the  run  length  Table  A2  (2  and  17).  The  addition  of  I at  0/70 
produces  agreement  in  both  directions  Therefore  go  to  the  bottom  of  the  next  column  (5)  and 
move  upwards  checking  the  row  totals.  At  40  5 the  addition  of  2 satisfies  the  two  totals.  In  the 
third  column,  however,  the  addition  of  I at  20  10  satisfies  column  10  but  still  leaves  a deficiency 
in  row  20  Since  the  following  flight  sequences  must  differ,  the  leading  diagonal  is  fixed  at  zero. 

Hence  row  20  must  wait  for  the  addition  of  I at  20  30.  The  process  continues  to  the  upper  right 
element  0 70.  with  all  totals  now  satisfied. 

The  resulting  table  is  only  one  of  many  possible  outcomes.  By  rearranging  the  ordering  of 
the  columns  and  rows  or  starting  at  a different  corner  of  the  table,  many  other  outcomes  would  . * 

result  Any  one  of  these  could  be  used  in  later  analysis. 


TABLE  A I 

Scaled  Flight  Sequence  Matrix 


TOF 

0 

5 

Following  TOF 

10  20  30 

40 

70 

Total 

0 

0 

10 

2 

1 

2 

0 

15 

5 

10 

9 

5 

1 

12 

1 

38 

10 

2 

9 

0 

2 

1 

0 

13 

20 

1 

5 

1 

0 

1 

0 

8 

30 

0 

I 

I 

0 

0 

2 

40 

3 

12 

1 

I 

0 

0 

0 

17 

0 

1 

0 

0 

0 

0 

1 

Total 

16 

38 

14 

8 

1 

16 

1 

TABLE  A2 

Scaled  Run  Length  Matrix 


i 


TOF 

1 

2 

3 

4 

5 

Run  Length 

6 7 8 

9 

10 

II 

12 

13 

Total 

0 

12 

3 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

17 

5 

16 

8 

7 

2 

2 

1 

2 

0 

1 

0 

0 

0 

1 

40 

10 

11 

4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

15 

20 

9 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

10 

30 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

40 

13 

4 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

19 

70 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

APPENDIX  2 
Matching  of  Spectra 


Suppose  we  have  a type  of  flying  (e.g.  TOFI)  with  a spectrum  as  shown. 

Exceedance  Levels  (g) 

-1-5  - 0-5  2 0 3 0 

k = [1  2 4 2] 

Then  if  N 100  flights  of  TOFI  the  spectrum  will  be: 

A k [100  200  400  200] 

Given  5 flight  trials  histories,  each  having  a spectrum  shown  below: 

Exceedance  Levels  (g) 

-1-5  -0-5  2 0 3 0 

A T2  I 2 1 

B 3 5 10  7 

C 0 14  2 

D 0 3 4 2 

E [5  10  10  5 

how  many  of  each  of  the  above  5 flights  will  produce  the  same  spectrum  as  N k? 

Let  Af.4  = No.  of  flights  of  type  A etc. 
with  x = [Xa.  Xu,  A 'c,  Xo.  A'e]. 

To  maintain  the  counts  at  each  level  we  get: 

"2  3 0 0 5 rAfal  p00‘ 

1 5 1 3 10  XK  200 

2 10  4 4 10  Xc  = 400 

I 7 2 2 5j  XD  l_200_ 

lxE_ 

or  briefly 

kft  x1  = N k' 

Also  we  require  a total  of  N 100  flights 
so  Xa  + Xu  |-  Xc  + Xt>  T Xk  = 100 

We  have  5 unknowns  and  5 equations,  so  our  first  choice  would  normally  be  to  use  any 
standard  method  for  solving  simultaneous  equations.  Fiowever,  we  have  a constraint,  namely 
that  no  value  of  X may  be  negative.  For  a best  solution  of  the  equations,  given  the  constraint, 
we  use  linear  programing. 

This  gives  the  rounded  results: 

x = [35  10  32  23  0], 

with  a total  of  100  flights  and  a spectrum  k„  with  contributions: 


L. 


Exceedance  Levels  (g) 
-1  -5  -0  5 2 0 3 0 


[diag  x]  kft  == 


No.  of 
Flights 


70 

30 

0 

0 


35  70 

50  100 
32  128 
69  92 


35 

70 

64 

46 


'35' 

10 

32 

23 

Totals  NA0  = k(nx  = [100  186  390  215]  I00 

- ,.Th'S  ,S  ,the  °p,l"1um  solutlon  Thc  discrepancies  present  can  be  corrected  by  altering  some 
of  thc  values  in  thc  actual  flight  histories. 

desert"'  * 'S  ‘I"  S°!Uri0n  °f  ,hC  m°re  gCneral  'inear  Pro8rarnming  problem  we  now 
describe.  The  spectrum  achieved  for  a given  trials  matrix  is: 

kft  * = Nk0,  say,  XA  + XB + XB  = N 

To  best  fit  the  required  spectrum  k for  the  given  TOF,  the  objective  function  to  minimize  is: 

J = Na‘|k  - k0|  = N*‘\kf,  x‘  - k| 

Where  a is  some  weight  vector  and  modulii  are  interpreted  termwise. 

The  constraints  of  course  are: 


x Ss  0 

and  XA  4 1-  Xg  = N 

In  the  example  a = [1  I 1 l )]. 


Computer  programs 

Program 

FATMET.F4 

DATRUN 

FTSEQ.F4 

MACEXT.F4 

TURNPT.F4 

FMSPEC.F4 

SIMPL 

MATAPE.F4 

MASEQ.F4 

TAPWRT 

FMKNT.F4 


APPENDIX  3 
List  of  Programs 

used  for  this  project  can  be  located  on  DECtape  358  and  accessed  using 
DPIP 


Function 

Produces  the  sequence  and  run  tables  for  the  raw  FM  data. 
Subroutines  called : SEQ 
Output  name:  PJHFM.OUT 

Scales  and  modifies  sequence  and  run  tables. 

Subroutines  called : ROUNOF,  FUDGE 
Output  name:  MODAT 

Produces  sequence  of  200  flights,  requires  program  RERUN.EXE 
on  disk. 

Output  name:  FLIGHT. SEQ 

Produces  50  channel  edited  data  from  flight  trials  data. 

Output:  To  magnetic  tape. 

Extracts  turning  points  of  Nt  from  edited  flight  trials  tapes. 

Output  name:  TURNPT.OUT 

Does  a fatigue  meter  type  count  and  simulates  quasi  flights  from 
TURNPT.OUT 

Output  name:  typed  in  manually. 

Matches  corrected  combinations  of  flight  trials  FM  counts  to  Pearce 
FM  data.  Uses  simplex  method  of  linear  programming. 

Output:  to  teletype. 

Produces  new  flight  sequence  using  results  of  SIMPL  with  original 
flight  sequence. 

Output  name:  FLTSEQ.NEW 

Fits  load  histories  to  each  of  the  200  flights  in  the  sequence. 

Output  name:  F00I  to  F200 

Transfers  200  files  in  order  from  disk  to  magnetic  tape. 

Does  a fatigue  meter  type  count  on  the  200  flight  histories. 

Output  name:  FM.OUT 
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